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Abstract 

We investigate quasi-periodic coronal pulsations during the decay phase of an X 
3.2 class flare on 14 May 2013, using soft X-ray data from the RHESSI satellite. 
Periodogram analyses of soft X-ray light curves show that ^ 53 s and ^ 72 s 
periods co-exist in the 3-6, 6-12 and 12-25 KeV energy bands. Considering the 
typical length of the flaring loop system and observed periodicities, we find that 
they are associated with multiple (first two harmonics) of fast magnetoacoustic 
sausage waves. The phase relationship of soft X-ray emissions in different energy 
bands using cross-correlation technique show that these modes are standing 
in nature as we do not find the phase lag. Considering the period ratio, we 
diagnose the local plasma conditions of the flaring region by invoking MHD 
seismology. The period ratio P 1 / 2 P 2 is found to be ^ 0.65, which indicates 
that such oscillations are most likely excited in longitudinal density stratified 
loops. 
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1. Introduction 


Solar flares are sudden release of magnetic energy occurring in the solar 
atmosphere, and lasting from a few tens of seconds (impulsive) to a few tens 
of minutes (gradual). Flare emissions are detected in the entire electromag¬ 
netic spectrum, ranging from radio, microwave, visible, ultraviolet. X-rays, 
hard X-rays, and even to Gamma-rays. The effect of these flares are also 
detected in the variations of the flux density of solar energetic particles ap¬ 
proaching towards the Earth. The electromagnetic radiation generated in solar 
flares often exhibits oscillatory patterns in the light curves. They are periodic 


....... ...... v^... 

ranging from a few milliseconds to several minutes ( 

Nakariakov & Verwichte. 

2OO5I: Nakariakov & Melnikov. 20091: Kim et al.l.l2012l: 

Kuorivanova et al.. 2013: 

Huang et al.l. 2014 and references cited thereinb These oscillations are referred 


to as quasi-periodic pulsations (QPPs). In some case s, they have inharmonic 
shape with apparent amplitude and period modulation (|Nakariakov fc Melnikov . 
20091) . 


QPPs are classified mainly as short QPPs (sub-second), medium QPPs (sec¬ 
onds to several minutes), and long QPPs (from several minutes to tens of 
minutes) according to the possible physical mechanisms associated wi th them 
( Nakyiakov fc Melnikovl . [2QQ9l:lKuDrivanova et al.LlioiQl) . According to lAschwandenI 
( 19871 ). the short QPPs are probably connected with the interaction of electro¬ 
magnetic, plasma, or whistler waves with accelerated particles detected in the 
radio emission. Medium QPPs are likely associated with magnetohydro dynamic 
(MHD) processes in the sol ar flaring loops as detected in radio, microwave, white 


light, and X-ray emissions (linglis 


Asai et aP. 2001 


et al 


ps as detected m radio, microwave, 
.120081 2009: Jakimiec fc Tomczakl. 




l2005:lReznikova et al.l. 12007: iMcAteer et al 


_ __iMelnikov et ah., . _ _ _ _ 

20051 : Reznikova fc Shibasaki , 201l[ and references cited therein). On the other 


hand, the l ong ones are usually relevant to act ive region (AR) dynamics (ISvch et al 


2009 . 2OI5I) and global oscillations of the Sun ( Tan et al. . l2010l:[Srivastava fc Dwivedil . 


2OIOD. 


QPPs in solar flares are one of the major diagnostic tools to study the 
physical conditions in flaring sites, and triggering mechanisms. These studies 
can also be extended to ste llar flares, which exh i bit QPPs in their radio, optical. 


and soft X-ray emiss i ons (Mathioudakis et al.l. 


IPandev fc Srivastaval . 12009^ Anflnogentov et al 


20031: Mitra-Kraev et afl 2005 . 

2013 : Srivastava et al.l . 2013[ ). 


Understanding QPPs during the flare eruption is still an open question. The 
fundamental physical processes of these oscillations in solar flares are relevant 
to understanding magnetic reconnection, magnetohydrodynamic (MHD) waves 
in coronal structures, particle acceleration, thermodynamics, and other kinetic 
effects. Although the properties of QPPs in solar hard X-ray (HXR) emissions 
have been widely studied by different researchers, their nature in solar soft X- 
ray emissions (SXRs) in different energy bands simultaneously remains to be 
extensively investigated. This calls for more observational information from 
samples with oscillatory behavior in low-energy X-rays. 

The ratio of the statistically significant QPPs of the simultaneously existing 
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spectral components is important for determining the nature and mechanisms of 


the pulsations both in 

observations and theorv (Van Doorsselaere et ah. 

2007; 

Srivastava et ah. 2008: 

Inslis & Nakariakov. 

20091 

Srivastava & Dwivedi. 

2010: 

Macnamara & Roberts 

. I 2 OII: Orza & Bahai 

, 2013 

: Luna-Gardozo et ah. 

2012; 

Erdelvi et al.L 2014: Guo et ah. 120151 and references cited there!. Different peri- 


ods of QPPs may be associated either with different MHD modes simultaneously 
present in a flaring; loop or with different harmonics of the same MHD mode 


present m a narmg; loop or witn dinerent narmomcs oi tUe same ivmjj mode _ 

(iMelnikov et ah . 2005 : Van Doorsselaere et ah . 2QQ7t Sriyastaya et alll2QQ8l:[Andries et al 


2QQ9I: linglis fc Nakariakovi I2QQ9I: ISrivastava fc Dwivedi |2Q1QI e.g.j. 

In this paper, we use the RHESSI data to study the intensity of the co¬ 
existing QPPs at 3-6 KeV, 6-12 KeV and 12-25 KeV multiple energy bands 
during the decay phase of an X-class solar flare on 14 May 2013. In Sect. 2, we 
briefly present the observations and power spectral analysis methods for detect¬ 
ing QPPs. Sect. 3 describes the results of the QPPs observed in different energy 
bands of the flaring loop. Exploration of various harmonics of appropriate MHD 
modes leaving their imprints on selected SXR while excited in the flaring loop 
system is discussed in Sect. 4. In Sect. 5, MHD seismology of flaring plasma 
using the observed QPPs are discussed along with conclusions. 


2. Observations and Analysis Methods 

The X 3.2 class solar flare on 14 May 2013 is investigated in the present 
study. It started at 00:42:39 UT, reached its maximum at 01:07:41 UT, gradually 
decayed and ended at 03:14:40 UT in GOES SXR observations. Its coordinate at 
the Sun is N11E74 in an active region NOAA 11748. The flare region is detected 
by both RHESSI and Nobeyama Radioheliograph (NoRH) in an integrated mode 
providing flaring light curves in its decay phase in X-rays and radio respectively. 
Eigures 1 (a,b,c) represent some features of this solar flare and related emissions. 

Eigure 1 (a) shows time profiles of the flare in soft X-rays taken with GOES satel¬ 
lite 1.0-8.0 A channel and microwave emission during this flare at 17 and 34 GHz 
observed by NoRH. NoRH provides intensity (Stokes parameter I=R+L) and 
circular polarization (Stokes parameter V=R-L) images at 17 GHz with a tem¬ 
poral resolution of 1.0 sec (Nakajima et al. 1994). The fluxes (Eigure 1(a), 
middle panel) and correlation curves (Eigure 1(a), bottom panel) at 17 GHz 
and 34GHz obtained from NoRH vary with the same pattern during the flare 
process. However, for the flare peak, the 17GHz flux is larger than the 34GHz 
and afterwards both fluxes become comparable with each other. The radio emis¬ 
sion is probably generated by gyro-synchrotron motion of accelerated electrons 
during the flare peak phase, while in the long decay phase the thermal freefree 
emissions dominate. The full Sun image obtained from NoRH at 17 GHz (Stokes 
parameter I=R-1-L) and the corresponding flaring region of the day under study 
is shown in Eigure 1 (b). 

Eigure 1(c)shows this localized temporal span of the light curve during the 
decay phase of the flare detected by RHESSI. In these data sets of 3-6 KeV, 6-12 
KeV and 12-25 KeV energy bands, best-fit exponential functions of the form / = 

Jo X are fitted for the time period between 01:23:40 UT and 01:51:56 UT 
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to remove the long-term background flare variations and corresponding graphs 
are shown in Figure 2.For 3-6 KeV energy band, the values Iq and b are 349.9 
and -0.003996; for 6-12 KeV data set, these values are 4265 and -0.004179; and 
for 12-25 KeV energy band, the corresponding values are 1960 and -0.005061 
respectively. 

We have studied the QPPs present in the de-trended signal during the decay 
phase of the flare by conventional Fast Fourier Technique (FFT), Maximum 
Entropy Method (MEM) and Lomb-Scargle (LS) periodogram analysis method. 
The MEM is an alternative method to EET which avoids the limited resolution 
and power 'leaking', due to the windowing of data, present in the EET, and it 
belongs to the class of methods which fit a satisfied model to the data. The 
parameters of a maximum entropy spectral estimation are equivalent to the 
' ones'’ in th e auto-regressive (AR) model of a random process in real domain 

( Biird.ll972h . _ _ _ 

On the other hand, the Scargle periodogram ( ScargleLll982l : lHorne fc Baliunad . 
is an important algorithm for time series analysis of unevenly sampled 


data. Being quite powerful for finding and testing, the significance of weak pe¬ 
riodic signals through false alarm-probability (EAP), a simple estimate of the 
significance of the height of a peak in power spectrum can be derived. In this 
m ethod, the confiden c e leve ls have been calculated according to the recipe given 
in iHorne fc Baliunasl ( 19861) . In case of EET and MEM spectral decomposition 


techniques, we have calculated the mean (/i) and standard deviation (cr) of the 
data points. We have then calculated the confidence levels as /r + 2cr = 95.4 
%; /i + 3cr = 99.7 % and /i -t- 4cr = 99.99 % respectively. X-axis of EET, MEM 
and LS graphs are drawn on a log scale. 

In the present study, the periods which are common in all the methods 
having power > 95 % confidence levels are considered only for final discussion 
and calculation. In order to justify the validity of the peaks found from the 
EET, MEM and LS methods, the confidence limits for different peaks o btained 
from the power spectra of different energy bands have b een calculated (|Haberl . 
19691: Chowdhurv fc Rav . 2006 : Chowdhurv et ah . 20091) . In every case, confi¬ 
dence levels are above 99%. In effect, we attempt to determine the interval in 
which any hypothesis concerning the periodicity of a certain solar event might 
be considered tenable and outside which any hypothesis would be considered 
untenable. The confidence limits (CL) are evaluated by generating a sample 
of 100 data points equally on both sides of a particular peak. This method is 
repeated for all the peaks under consideration within a spectrum. The peak 
that is sharp gives the minimum value of standard error and closer values of 
confidence limits. 

The confidence interval provides the lower and upper limits to which the 
population parameter has a high probability of being included. The population 
parameter standard deviation a can be calculated from the following formula: 

a={ETi/{N-l)-{ET^/N{N-l)f'^ ( 1 ) 

Here Ti represents the i value of the data under study and N is the total 
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number of data points. The standard error SE^ is the standard deviation of 
the sample mean (from sampling distribution) which is estimated as: 

SEm = a/{Nf^ ( 2 ) 

The confidence limits (CL) for 99% confidence can be computed as 

CL = T,,e±2.58(5E;^), (3) 

where Tave is the mean value of the time period of the sample data points. 

To determine confidence limits we have used the formula of Za/2 x (n). 
Here Za/2 is the confidence coefficient, a is confidence level, a is standard de¬ 
viation, and n is sample size. Then, we have converted the 99 % to a decimal 
0.99 and divide it by 2 to get 0.495. According to the table the closest value 
corresponding to 0.495 is 2.60. Then we have multiplied 2.6 by 0.99 (critical 
value by our standard error) and we get 2.574 (^2.58) and this is our margin 
of error. 

We have then calculated the cross-correlation coefficients (CC) between de¬ 
trended data of three different energy channels observed by RHESSI to study the 
time-delay with various time-lags (0, ±4, ±8...± 50 seconds). Each of these data 
sets has been shifted forward or backward and the cross-correlation coefficient 
has been calculated after shifting the time-lag (L) corresponding to the best 
correlation of the time phase under investigation. The results are shown in 
Eigures 6 (a,b,c). 

3. Results 

3.1. Periodogmm Analysis 

Eigure 3(a, b, c) represent the EET, MEM and Scargle spectrum of 3-6 KeV 
RHESSI channel and the corresponding results are given in Table 1. 

Eigures 4 (a, b, c) represent the EET, MEM and Scargle spectrum of 6-12 
KeV RHESSI data and the corresponding results are given in Table 2. 

Eigures 5 (a, b, c) represent the EET, MEM and Scargle spectrum of 12-25 
KeV RHESI data and the corresponding results are given in Table 3. Thus 
we find, in all methods and in all energy bands, the common periods to be 
^71 sec (Pi) and ^53 sec (P 2 )- However, for the energy band of 12-25 KeV, an 
additional period ^39 sec (P 3 ) is found in all spectral decomposition techniques. 
The period ratios are estimated as P 1 /P 2 = 71/ 53 ^ 1.339; Pi/2P2^0.669; 
Pi/P 3=71/39-1.82 or P 3 /P 1 -O. 549 . 

3.2. Cross-correlation Analysis 

We also study lagging and leading behaviors of RHESSI 3-6 KeV, 6-12 
KeV and 12-25 KeV channels during the decay phase of the flare by the cross¬ 
correlation technique and the plots are displayed in Eigures 6 (a, b, c). It is 
detected that there is no lag/lead between different energy channels demonstrat¬ 
ing very high correlation of the signals in all observational bands. Erom figure 
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Table 1: List of periods detected with their Standard Errors for Confidence Limits in the 
decay phase of solar fiare in energy band 3-6 KeV 


Energy band 


Identified spectral peaks (sec) 



1 

2 

3 


(a) By FFT 


Tave 

71.56 

53.34 

38.62 

SE„, 

0.035 

0.0195 

0.0101 

T 

71.56 ± 0.0903 

53.34 ± 0.0503 

38.62 ± 0.0261 


(b) By MEM 


Tave 

53.36 

71.74 

32.71 

SE„, 

0.274 

0.457 

0.102 

T 

53.36 ± 0.7069 

71.74 ± 1.1791 

32.71 ± 0.2632 


(c) By Scargle Technique 


Tave 

71.89 

53.47 

38.80 

SE„, 

0.0352 

0.0192 

0.0101 

T 

71.89 ± 0.0908 

53.47 ± 0.0495 

38.80 ± 0.0260 


Table 2: List of periods detected with their Standard Errors for Confidence Limits in the 
decay phase of solar fiare in energy band 6-12 KeV 


Energy band 


Identified spectral peaks (sec) 

1 

2 

6-12 KeV 

Tave 

(a) By FFT 

71.96 

53.34 


SE„, 

0.0363 

0.0203 


T 

71.96 ± 0.0937 

53.34 ± 0.0524 


Tave 

(b) By MEM 

71.80 

53.43 


SE„, 

0.479 

0.2803 


T 

71.80 ± 1.2358 

53.43 ± 0.7232 


Tave 

(c) By Scargle Technique 
71.73 

53.41 


SE„ 

0.0354 

0.0183 


T 

71.73 ± 0.0913 

53.41 ± 0.0472 
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Table 3: List of periods detected with their Standard Errors for Confidence Limits in the 
decay phase of solar flare in energy band 12-25 KeV 


Energy band 


1 

Identified spectral peaks (sec 
2 

) 

3 

4 

12-25 KeV 

Tave 

71.08 

(a) By FFT 

53.14 

38.88 



SE„, 

T 

Tave 

0.0312 

71.08 ± 0.0903 

71.52 

0.0185 

53.34 ± 0.0503 
(b) By MEM 

53.36 

0.0106 

38.88 ± 0.0273 

32.69 

39.76 


SE„, 

0.537 

0.278 

0.104 

0.213 


T 

71.52 ± 1.385 

53.36 ± 0.7172 

32.69 ± 0.2683 

39.76 ± 0.5495 


Tave 

SE„, 

T 

72.15 

0.0367 

72.15 ± 0.0947 

(c) By Scargle Technique 

53.60 

0.018 

53.60 ± 0.0464 

38.78 

0.0104 

38.78 ± 0.0268 



6 (a), it is detected that maximum cross-correlation coefficient is 0.68±0.0168 
at 0 second. Figure 6 (b) exhibits that maximum cross-correlation coefficient 
is 0.643±0.0195 at 0 second. The maximum cross-correlation coefficient is 
0.781±0.0264 at 0 second between the energy range 6-12 KeV and 12-25 KeV 
(Fig. 6c). 

3.3. Morphology of the Flaring Loop System 

GOES SXI image of the active region and its loop system beneath where the 
flare occurred is shown in Fig. 7 (right panel), while the same active region in 
17 GHz NoRH observations is also displayed in Fig. 7 (left panel). The North- 
South axis in the GOES map is inclined for aproximately 20 degrees clockwise 
relatively to the vertical direction. Therefore NoRH map shows the true position 
of the AR while GOES SXI image shows the rolled partial view. The soft X- 
ray enhancement occurred in the loop system during the evolution of the solar 
flare. The length of the flaring loop system is derived as L^190 Mm while its 
width is ^12 Mm. The magnetic configuration of the flaring site obtained for 
01:00 UT when approximate loop-length and width are derived, remains almost 
the same and changes insignificantly in 20 min later when QPOs evolve in the 
X-ray light curves. The loop-length is estimated by tracking the pixels along 
the loop. However, it should be noted that the estimated length is the lower 
bound projected length. The loop width is also estimated across the loop near 
its apex. 


4. Results and Discussion 

The observed periods are most likely associated with the first two harmonics 
of the fast magnetoacoustic wave modes in the flaring loop system. The phase 
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speed of the fundamental mode of magnetoacoustic wave is Vph = 2LjP^ 5200 
km/s, which seems to he with the fastest magnetoacoustic branch (tubular¬ 
mode) within the loop. Since the QPPs modulate the emissions coming from 
the flaring loop system, they are compressible fast modes known as fast magne¬ 
toacoustic sausage modes. 

The dispersion relation of fast-sausage waves in straight cylindrical tube 
includes both tra pped and leaky modes de pending on the plasma and magnetic 
field parameters ( Edwin fc Roberts . 19831) . 

The cut-off wavenumber kc is written as given bv lEdwin fc RobertsI ( 1983 ) 


and [Roberts et al.l (|l984l) 


k = kc 


VI 


VI 


VI 


1/2 . 
Jo 


( 4 ) 


where Ya and Yac are Alfven speeds inside and outside the loop, a is the 
loop radius and jo= 2.4 is the first zero of the Bessel function. The modes with 
k > kc are trapped in the loop system, while the modes with k < kc are leaky. 
The trapped fundamental global sausage mode (P = 72 s) can be realized in 
this loop system only if the density contrast is: 


II < 

2jo Eao 


0.65, 


( 5 ) 


which is equal to 593 for a given loop length L^190 Mm and width a^6 
Mm. Therefore, the observed flaring loop must be very thick loop with sig¬ 
nificant density contrast to possess the global sausage modes. The maximum 
phase speed of the sausage wave 5200 km/s is considered as an external Alfven 
speed (VyLe), while Yao is the internal Alfven speed within the denser loop sys¬ 
tem. Using Eq. (2), we can estimate the ratio of external and internal Alfven 
speed (V^e/VAo) in the loop system as ^24.0. The phase speed of the funda¬ 
mental mode sausage wave {Yph) can be maximum limit to the external Alfven 
speed, i.e., VAe=Vp/i^5200 km/s. Therefore, the internal Alfven speed {Yao) 
within the loop will be 216 km/s. The diagnosed internal Alfven speed {Yao) 
of the flaring loop system is 216 km/s. The Alfven speed inside the bulky flar¬ 
ing loop is significantly less compared to that in the outside ambient medium. 
This is expected as the density of such loops is higher compared to the normal 
loops as well as ambient plasma. This loop system does not possess the global 
sausage mode s with period greater than PGSM^27r/kc Yac^ 2L/VAe^73.0 s 


(Pascoeetal 


2QQ7I ). It should be noted that L=7r/ kc. The longest period of 
the trapped sausage mode can have the period of ^73 s. The observed period of 
fundamental mode sausage waves is less than this maximum period. Therefore, 
this period and its harmonics are well accepted in the given morphology and 
density contrast of the loop system. 

Eurther, we have also found zero time lag between different soft X-ray en¬ 
ergy bands. This zero lag in all bands depicts the periodicity that the wave 
concurred by sausage mode is not a propagating wave, but it is standing wave 
in nature. Therefore, we basically observe fast MHD pulsation of short periods 
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excited due to flare activity in the flaring loop system. As a result, the observed 
denser loops can support the non-leaky, trapped glob al sausage mode in the 
estimated length and width as well as density co ntrast ( Nakariakov et all . I2QQ3I: 
lAschwanden et ahl . |2QQ4 ISrivastava et al.l . 2008). Such flaring/post-flare loops 
usually have a very high density contrast of the order of 10^-10^. The density 
evolves very rapidly during the flare heating and makes such loops highly denser. 
The observed loops are the X-class flaring loops where large heating is likely to 
take place during the strong X-class flare. Therefore, the estimated high den¬ 
sity contrast is likely to occur and also falls in the expected range. Thus, the 
trapped sausage mode is probably excited in the loop system. 

The ratio between the periods of the fundamental and the first harmonics of 
sausage waves is equal to Pi/P2^1.34, which is shifted from 2.0. Similar phys¬ 


ical behavior is observed for the fast-kink oscillations (IVerwichte et al.L I2QQ4 ; 


Van Doorsselaere et ah . ^QQ^jl: IPe Moortel fc Bradv . 2007: IWhite fc Verwichte . 


2 OI 2 I: Srivastava fc Goossensl. 20131). for sausage waves ( Srivastava et al.l . 2008 ) 
and for slow waves ([Srivastava fc Dwivedi , 2010[) . The deviation of Pi /P 2 from 


2.0 in homogeneous loops i s very small due to the wave dispersion (|McEwan et al 


20061: lAndries et al.l . 12005 ). Howeve r, longitudinal density stratification causes 
significant shift of P 1 /P 2 from 2.0 (Andries et ah. 2005. 20091: McEwan et al 


2006: Van Doorsselaere et al. ll2nn7l:M 


acnamara . 


fc Roberts. 2010[|2011: Lnna-Cardozo et al 


2 OI 2 I and references therein). The longitudinal variation of density in a coronal 
loop is the most likely cause for the departure of the observed period ratio of 
the sausage waves. How ever, magnetic field geometry may also play an impor¬ 
tant ro le in its deviation ( Macnamara fc Robertsl . [ 2 OII ). iMacnamara fc RobertsI 
(|201l[ ) have shown in their analytical calculations that as V^e/VAo tends to 
infinity the period ratio of the sausage waves (P 1 / 2 P 2 ) approaching towards 
0.5 in their slab-model. Normally, it should be equal to 1.0. Taking account 
of our observed parameters, we obtain a very high (^24) external to internal 
Alfven speed ratio. The period ratio (P 1 / 2 P 2 ) is significantly less, i.e., 0.65. 
This observation is qualitativ ely in agreement with the analytical solution of 
Macnamara fc Robertsl ( 201l[ ) that period ratio reduces significantly for higher 
Alfven velocity ratios. 

In the present work, we find that QPPs are most likely caused by the sausage 
modes, and therefore, based on their signature local plasma properties are di¬ 
agnosed. However, there may be other possible mechanisms which may be 
potentially responsible for QPP modulation. A few examples may be cited as 
the periodic variations of the angle b etween the line-of-sig ht and magnetic field 
by weakly compressible kink waves ( Cooper et al.l. 2QQ3I). or by excitation o f 
QPPs by MHD oscillation of a n external loop (|Nakariakov fc Melnikovl . l2QQ9[ ). 
Recently. iKolotkov et al.l ( 2015 ) have analyzed the multiple MHD modes during 
the impulsive phase of the same flare that was approximately 18 min before the 
time interval analyzed in our present paper and lasting for a few minutes only. 
The signatures of waves are revealed with the mean periods of 15, 45, and 100 
s. The 15 s and 100 s period are associated with the MHD modes rapidly de¬ 
caying within a few wave periods. They are interpreted as leaky sausage modes 
and the kink modes respectively. The 45 second periodicity showed rather non 
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decaying behavior and the period similar to the findings of our present paper. 
This example shows that in the flaring region, the plasma and magnetic field 
conditions are rapidly varying. Therefore, even a single flare energy release over 
different epoch in the given duration of the flare may give rise to the evolution of 
different MHD modes and their properties. These properties are highly linked 
to the local plasma and magnetic field conditions, and can be used in diagnosing 
highly variable atmosphere there. 
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May 14, 2013 / 02:44 UT / NoRH, R+L 



Figure 1: (a) GOES plot showing temporal evolution of soft X-rays (top panel). The mi¬ 
crowave emission measured by NoRH at 17 and 34 GHz - integral flux (middle panel) and 
correlation plots (bottom panel), (b) Full Sun in intensity NoRH channel at 17 GHz. The 
flare region is indicated by the solid box. (c) Localized temporal span of the light curves of 
different energy bands during the decay phase of the ffare detected by RHESSI on May 14, 
2013 
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Figure 2: RHESSI soft X-ray light curves on May 14, 2013, from 01:23:40 to 01:51:56 UT. (a) 
Blue line is the trend htting X-ray light curve for 3-6 Kev data; (b) same for energy band 
6-12 keV and (c) same for energy band 12-25 Kev. 



Figure 3: (a) FFT spectra of RHESSI 3-6 Kev data; (b) MEM spectra for the same data and 
(c) Scargle spectra of the same data. 



Figure 4: (a) FFT spectra of RHESSI 6-12 Kev data; (b) MEM spectra for the same data 
and (c) Scargle spectra of the same data. 
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Figure 5: (a) FFT spectra of RHESSI 12-25 Kev data; (b) MEM spectra for the same data 
and (c) Scargle spectra of the same data. 



Figure 6: (a) Cross-correlation analysis between RHESSI 3 -6 and 6-12 Kev data sets ; (b) 
Cross-correlation analysis between RHESSI 3-6 and 12-25 Kev data ; (c) Cross-correlation 
analysis between RHESSI 6-12 and 12-25 Kev data set. 
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Figure 7: Left-panel : NoRH image showing the position of flaring active region. Right- 
panel : Partial view of the GOES SXI image of X-class flaring region (negative color image) 
and loop system beneath which the flare occurred. The soft X-ray enhancement occurred 
in the loop system is shown during the evolution of the solar flare. The North-South axis 
in the GOES map is inclined for aproximately 20 degrees clockwise relatively to the vertical 
direction. Therefore NoRH map shows the true position of the AR while GOES SXI image 
shows the rolled partial view. 
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